Ultra-fast quenching of binary colloidal suspensions in an external magnetic field 
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An ultra-fast quench is applied to binary mixtures of superparamagnetic colloidal particles con- 
fined at a two-dimensional water- air interface by a sudden increase of an external magnetic field. 
This quench realizes a virtually instantaneous cooling which is impossible in molecular systems. 
Using real-space experiments, the relaxation behavior after the quench is explored. Local crystal- 
lites with triangular and square symmetry are formed on different time scales and the correlation 
peak amplitude of the small particles evolves nonmonotonicaliy in time in agreement with Brownian 
dynamics computer simulations. 
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Temperature quenching belongs to the key processing 
techniques to produce amorphous and crystalline solids 
which are considerably different from their thermody- 
namically stable counterparts. Possible applications can 
be found in metallurgy, ceramics and semiconductor dop- 
ing. For example, quenching is used to construct ceramic 
material with a high mechanical stability ll| , to design 
resistance devices for integrated circuits [2j, and to im- 
prove the optoelectronic properties of semiconductors . 

Temperature quenching techniques are most efficient if 
they are performed suddenly, i.e. if the system tempera- 
ture changes on a time scale that is much shorter than a 
time upon which a typical particle motion is performed. 
While this can be realized in computer simulations, see 
e.g. [3, 0], it is practically impossible to be achieved 
for molecular systems where the quench is performed by 
a coupling to an external heat bath. There, it takes 
some time until the prescribed temperature is realized 
in the sample. However, as we shall show in this Letter, 
an ultra-fast quenching is possible for colloidal particles 
which move much slower and are highly susceptible to 
external fields. 

If a suspension of superparamagnetic particles is used, 
an external magnetic field induces magnetic dipoles in the 
particles which gives rise to dipolar interactions between 
them [7] . In equilibrium, temperature and magnetic field 
strength determine the dimensionless coupling strength 
between the particles such that temperature is strictly 
equivalent to the inverse square of the magnetic field 
strength @, Therefore a sudden increase of the ex- 
ternal magnetic field corresponds to an ultra-fast quench 
towards lower temperature. It is important to note that 
the increase of the external magnetic field occurs on a 
time scale of a few ms much smaller than a couple of 
seconds typically needed by a colloidal particle to diffuse 
over its own size. Thereby - unlike molecular systems - 
colloidal systems can be exposed to ultra-fast tempera- 
ture quenches. 

In this Letter, we exploit this idea of quasi in- 



stantaneous quenching for a binary suspension of two- 
dimensional superparamagnetic colloidal particles con- 
fined to a planar water- air interface. At high temper- 
atures - or equivalently at low external magnetic field 
strengths - the system is weakly correlated. After an 
abrupt increase of the external magnetic field, the re- 
sponse of the system and the early stage relaxation to- 
wards its new state is monitored by real-space imag- 
ing. For the prescribed composition of the mixture, the 
equilibrium state is a crystalline lattice with alternat- 
ing stripes of pure triangles of the majority component 
and mixed squares [10[. This complicated stable crystal 
is not reached on the time scale explored but the sys- 
tem reveals structural heterogeneities corresponding to 
local metastable patches of crystalline order. These crys- 
talline zones are forming on different time scales which 
also gives rise to nonmonotonic behaviour in time for 
structural correlation peaks of the small particles. The 
real-space experiments are in agreement with our Brow- 
nian dynamics computer simulations. Our results can be 
used to steer the microstructure of composite materials 
upon quenching and reveal the interplay between vitrifi- 
cation and crystal nucleation [ll|. 

The experimental system consists of a suspension of 
two kinds of spherical super-paramagnetic colloidal par- 
ticles denoted by A and B. Those particles are char- 
acterized by different diameters (cIa = 4.5 /im, cIb = 
2.8 /im) and magnetic susceptibilities (xa = 6.2 x 
lO" 11 Am 2 /T,XB = 6.6 x 10 -12 Am 2 /T). The relative 
composition X = Nb/(Na + Nb) of B particles is fixed 
at 40%. Due to their high mass density, the particles are 
confined by gravity to a flat water-air interface formed 
by a pending water drop. The droplet is suspended by 
surface tension in a top sealed cylindrical hole (6 mm di- 
ameter, 1 mm depth) on a glass plate. The system can 
be considered as ideally two dimensional since the ther- 
mally activated "out of plane" motion of the particles is 
of the order of a few tens of nanometer. A coil produces a 
magnetic field H perpendicular to the water-air interface 
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inducing a magnetic moment (ui z = XiH with i = A, B) 
in each particle which leads to a repulsive dipole-dipole 
pair interaction of the form 

u ij (r) = ^ XiXj H 2 /r 3 (i,j = A,B). (1) 

For this inverse power potential, at fixed composition X 
and susceptibility ratio Xb/xa, all static quantities de- 
pend solely [lH on a dimensionless interaction strength 
(or coupling constant) 

4irk B Ta 3 1 } 

where k B T is the thermal energy at room temperature 
and a = 1/ yfpA the average interparticle separation be- 
tween A particles. The partial A particle density is set 
to pa = 1.97 x 10" 3 /im -2 so that a ~ 22.5/mi. The pair 
interaction is directly controlled over a wide range via the 
magnetic field. Making use of video microscopy, trajec- 
tories of all particles in the field of view can be recorded 
over several days providing the whole phase space infor- 
mation. The quench was realized upon suddenly rising 
the coupling T from 1 to 71 at time t = by increasing 
the magnetic field. The time scale of the quench is only 
limited by the electronics of power supply and was mea- 
sured to be 5 ms. This is much faster than the typical 
Brownian time r = (1 2 a /Da needed for an A particle to 
diffuse over its own distance with Da = 0.11/im 2 /s be- 
ing the short-time diffusion constant for A particles. The 
time r is measured to be 50 s. Hence the quench can be 
considered to be ultra-fast. In a unquenched system clear 
deviations of purely diffusive behavior appear already at 
T = 30. For T = 71 a well pronounced plateau between 
t/r = 2 and t/r = 200 is observed in the mean square 
displacement with an inflection point at t/r = 40 [7[. 
Hence the quench can be considered to be ultra-fast and 
deep into the supercooled stage. 

In parallel, we perform nonequilibrium Brownian dy- 
namics (BD) computer simulations [l3[ of our experimen- 
tal system described above. We employ point-like dipoles 
that interact following Eq. ([T]) with xb/xa — 10% and 
X = 40% in accordance to the experimental parameters. 
Knowing that the diffusion constant scales with the in- 
verse of the radius of a particle, D B was chosen such that 
D B /D A = d A /d B . N A = 400 A particles and N B = 267 
B particles were placed in a square box with periodic 
boundary conditions applied in the two directions. A 
finite time step St = 6 x 10 _4 r was used. The early 
stage of the system response to the ultra-fast quench has 
been observed both in real-space experiments and BD 
computer simulations and was characterized by various 
time-dependent correlations. Our first aim is to iden- 
tify the dynamical formation of local crystallites after 
the quench and detect building blocks of the underly- 
ing equilibrium crystal. The latter consists of alternat- 
ing stripes of pure A triangles and intersecting squares 




FIG. 1: (a) Fraction of B particles belonging to a crystalline 
square surrounding (see inset) and (b) fraction of A particles 
belonging to a crystalline triangular surrounding (see inset) 
versus reduced time t/r for an ultrafast quench from r = 1 to 
r = 71. The lines are experimental data while the symbols (*) 
are data from BD simulations. Two experimental snapshots 
for a time t/r — 0.6 just after the quench (left configuration) 
and a later time t/r = 60 (right configuration) are shown. 
Big particles are shown in blue if they belong to a triangular 
surrounding and in red if the belong to a square surrounding. 
All other big particles are shown in white color. Few big parti- 
cles belonging to both triangular and square surroundings are 
shown in pink color. The small particles are shown in green 
if they belong to a square center of big particles, otherwise 
they appear in yellow. Also included are BD simulation data 
for an instantaneous "steepest descent" quench from r = 1 
to r = oo (A) and for a linear increase of T from r = 1 to 
r = 71 on a time scale of 30r (+). 



of A and B particles [10]. Therefore we have used cri- 
teria to define A particles which have a pure triangular 
surrounding of other A particles, i.e. which are close to 
a cut-out of a pure triangular A crystal, and, likewise, 
we have identified A and B particles which form locally 
an equimolar square lattice S(AB)[14\. The correspond- 
ing two structure elements are shown as insets in Figure 
[TJ In detail, we associate a triangular surrounding to 
an A particle if the following two criteria are fulfilled 
simultaneously [§]: 1) The 6-fold bond order parame- 
ter p 6 = V^p^ (where ^ 6 = \ YInn ex P (^nn) with 
6nn denoting the angles of the six nearest neighbour 
bonds relative to a fixed reference) is larger than 0.94. 2) 
The relative bond length deviation b$ = - YInn 
where I is the average length of the six bond lengths Inn 
is smaller than 0.04. This double condition selects local 
configurations close to those of a perfect triangular lat- 
tice where p$ is unity and b§ vanishes. Likewise we define 
a square surrounding around a B particle by the criteria: 
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1) The 4- fold bond order parameter^ = ^^4^4 (where 
^4 = I XIatat ex P (^nn) with Onn denoting the bond 
angles of the four nearest neighbour AB bonds) is larger 
than 0.92. 2) The corresponding relative AB bond length 
deviation 64 is smaller than 0.05. 

Experimental snapshots before and after the instan- 
taneous quench are shown in Figure [T] with color-coded 
particles indicating the locations of A and B particles 
which belong to local triangular and square clusters. If 
an A particle has a triangular surrounding all seven A 
particles including the central one with its full surround- 
ing are shown in blue. Conversely, if a B particle has a 
square surrounding it is colored in green and its four A 
neighbours are colored in red. All particles which belong 
both to the blue and red class are shown in pink. From 
the snapshots of Fig. HJ it is evident that the crystalline 
clusters form locally and grow as a function of time. Pref- 
erentially triangular structures form in an area depleted 
from B particles while square structures nucleate around 
"seeds" which possess a structure close to an underly- 
ing square. The resulting crystalline patches are then 
fluctuating with a life-time of about 30r. The fraction 
of A particles with a triangular and of B particles with 
a square surrounding are shown as a function of time 
in Figure [T] a) and b), respectively. On the time scale 
considered in this plot one finds an increase from almost 
zero before the quench to 12 percent for the triangular 
A particles and the square B particles. The time scales 
upon which triangular and square structures are formed 
are different by a factor of about 2. The triangular struc- 
ture is forming more rapidly than the square one. This 
follows from the fact that an AB square structure re- 
quires more fine-tuning of structural correlations of both 
species than a triangular one which can directly emerge 
in regions depleted from small particles. The number 
of pink A particles which belong to both triangular and 
square surroundings is growing on a time scale slightly 
longer than that for the square structure to a fraction of 
only 1-2 percent far away from the equilibrium structure 
where the fraction of pink particles is 2/3. Last but not 
least, within the statistical error, there is agreement be- 
tween the experimental data and our nonequilibrium BD 
computer simulations. 

After a relaxation time of t = lOr, we have further 
calculated the mean-square displacements for particles 
which are crystalline (i.e colored in blue, red, green) and 
non-crystalline particles during a time of 60r, see Figure 
1. In both experiment and simulation, for both particles 
species the mean-square displacement of the uncolored 
particles is twice as large as that of the colored ones. 
This gives clear evidence for a correlation between local 
slow dynamics and local crystallinity. 

By simulation we have finally explored a different 
depth and rate of the quench (see Fig. 1). While a "steep- 
est descent" quench leads to a faster growth of crystalline 



patches and to an almost doubled crystalline fraction of 
particles, a linear increase of V within a time window of 
30r delays their formation accordingly. 
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FIG. 2: Partial pair distribution functions gAA(r),gAB(r) and 
gBB(r) of A- and B-particles versus reduced distance r/a at 
three different reduced times (a),(d),(g) t/r = 0; (b),(e),(h) 
t/ T = 1; (c),(f),(i) t/r = 80. BD results (solid lines) are 
compared to experimental data (symbols). 

Next we show the time evolution of the partial pair dis- 
tribution functions ^(r), 9ab{t) and ^(r) for three 
times t/r = 0, 1, 80 in Figure O First of all, there 
is again agreement between simulation and experiment 
While the correlations in the starting configurations be- 
fore the quench are weak, they are quickly increasing as a 
function of time towards a strongly correlated glass. An 
averaged square order can be extracted from the growth 
of an intermediate peak in #aa(0 at r = \[2a. This 
peak grows much slower than the first peak amplitude. 
The growth of the first peak amplitude in gAAir) and 
gAB(r) are monotonic in time, whereas there is a non- 
monotonicity in that of ^(r). This is clearly visible 
in Figure [3] where the dynamical evolution of the am- 
plitudes of all three partial pair distribution function is 
shown. With in the statistical accuracy, the peaks of 
9aa(^) and #ab(0 grow on the same time scale in a 
monotonic way, while the peak of qbb^) overshoots its 
final equilibrium limit, both in experiment and simula- 
tion. 22]. We explain this striking effect by a two-stage 
relaxation process of the small particles which are first 
excluded from the triangular crystallites regions of the 
big ones. Concomitantly they show a strong correlation 
since they are compressed until the total system relaxes 
back to a state where the small particle optimize their 
correlations in the network dictated by the big ones. 

The energetic optimization to the final state is shown 
in Figure [4] revealing agreement between experiment and 
simulation within the statistical inaccuracies. There is a 
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FIG. 3: Amplitude g m ax of the first correlation peak in the 
partial pair distribution function qaa {t) (triangles), gAB(r) 
(crosses), gBB(r) (stars) as a function of reduced time t/r. 
Note the non-monotonicity in the amplitude of gBB(r). (a) 
Brownian dynamics data, (b) experiments. 



huge jump in the potential energy per particle immedi- 
ately after the quench which then relaxes quickly towards 
a transient state. Then, on a second time scale, we find 
a slower relaxation process accompanied by a structural 
ordering as can be seen in the corresponding experimen- 
tal snapshots right before and after the instantaneous 
quench (see inset of fig. 



good agreement with Brownian dynamics computer sim- 
ulations. 

Our real-space characterization allows to identify the 
pathways of relaxation into a quenched glassy phase. As 
the quench is ultra- fast, the dynamical pathways are not 
blurred by an additional time-scale from the quench his- 
tory. The built-up and the fluctuations of local crystal- 
lites after the quench can be directly followed. There- 
fore our analysis can help in a more fundamental way to 
understand the interplay between vitrification and crys- 
tallization In fact, the structural heterogeneities 
detected here give a considerable weight to concomitant 
dynamical heterogeneities [TBI, [l6[ and therefore represent 
a relevant contribution to the dynamically heterogeneity 
of glasses [13, [H| . Finally, our results may also be useful 
to identify pathways of defect annealing in the crystalline 
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FIG. 4: Time evolution of the potential energy u per particle 
in units of UbT versus reduced time t/r. Brownian dynamics 
simulation data are shown in red, experimental data in green. 
For a better comparison the jump after the quench is trun- 
cated (see the cut in the y-axis). The inset is an expanded 
view. Corresponding experimental snapshots are shown be- 
fore the quench at t/r = 0~ on the left side and after the 
quench at t/r = 80 on the right side. Big particles are shown 
in white while small particles are presented as yellow spheres. 

In conclusion, we have realized an ultra-fast quench in 
two-dimensional colloidal mixtures from a weakly inter- 
acting to a strongly interacting situation by a sudden in- 
crease of an external magnetic field which controls the in- 
terparticle repulsion. The system spontaneously relaxes 
by exhibiting structural inhomogeneities which reflect the 
underlying stable crystal and correlate with slower re- 
gions in the dynamics. The experimental data are in 
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